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ABSTRACT 

We discuss a mechanism for cosmic ray penetration into an interplanetary 
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^ . magnetic flux rope, particularly the effect of the finite Larmor radius and laag- 

er . 

netic field irregularities. First, we derive analytical solutions for cosmic ray be- 
havior inside a magnetic flux rope, on the basis of the Newton-Lorentz equation 
of a particle, to investigate how cosmic rays penetrate magnetic flux ropes under 
an assumption of there being no scattering by small-scale magnetic field irregu- 
larities. The results show that the behavior of a particle is determined by only 



r^ , one parameter /o; that is, the ratio of the Larmor radius at the flux rope axis 

O I to the flux rope radius. The analytical solutions show that cosmic rays cannot 

penetrate into the inner region of a flux rope by only gyration and gradient- 
curvature drift in the case of small /q. Next, we perform a numerical simulation 
k> ! of a cosmic ray penetration into an interplanetary magnetic flux rope by adding 

small-scale magnetic field irregularities. The results show that cosmic rays can 
penetrate into a magnetic flux rope even in the case of small /o because of the 
effect of small-scale magnetic field irregularities. This simulation also shows that 
a cosmic ray density distribution is greatly different from that deduced from 
a guiding center approximation because of the effect of the finite Larmor radius 
and magnetic field irregularities for the case of a moderate to large Larmor radius 
compared to the flux rope radius. 

Subject headings: cosmic rays — diffusion — interplanetary medium — solar- 
terrestrial relations — Sun: coronal mass ejections (CMEs) — turbulence 
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Introduction 



Sudden decreases in cosmic ray intensity during geomagnetic storms were first observed 
by Forbush (1937) and are now well known as Forbush decreases. Observations of cosmic 
rays by neutron monitors and muon detectors revealed that Forbush decreases are caused by 
the passage of interplanetary shock waves with downstream highly turbulent magnetic fields 
and by magnetic clouds that have well-ordered strong magnetic fields. The observations 
also showed that there is an anisotropy in the cosmic ray distribution prior to and during a 
cosmic ray intensity decrease (Cane 2000 and references therein). 

When interpreting and modeling the observations, it is very important to clearly sepa- 
rate the shock wave effect and magnetic cloud effect because their characteristics are quite 
different (Wibberenz et al. 1998). There are many theoretical and numerical models for the 
shock wave effect (e.g., Nishida 1982; Thomas & Gall 1984; Chih & Lee 1986; le Roux & 
Potgieter 1991). Although data analyses of the magnetic cloud effect have also been carried 
out (e.g., Sanderson et al. 1990; Ananth & Venkatesan 1993; Cane et al. 1994; Ifedili 2004; 
Subramanian et al. 2009), only a few theoretical and numerical models of the effect have 
been studied (Cane et al. 1995; Munakata et al. 2006; Kuwabara et al. 2009), and this has 
limited our understanding of the behavior of cosmic rays inside magnetic clouds. 

Cane et al. (1995) (see also Wibberenz et al. 1998; Cane et al. 2000) modeled a cosmic 
ray penetration into a magnetic cloud by cross-field diffusion. They expressed the penetration 
by means of diffusion equation in a static vacuum cylinder region and deduced the cross-field 
diffusion coefficient for ^ 1 GV cosmic ray protons. Munakata et al. (2003) and Kuwabara 
et al. (2004) applied the muon observation data to deduce magnetic cloud geometry. These 
were the first derivations of the three-dimensional geometry of magnetic clouds from muon 
observation data. They assumed that the cosmic ray density distribution inside the magnetic 
cloud is Gaussian. Validity of the Gaussian distribution was not mentioned in the paper 
although the results were in good agreement with the geometry deduced from interplanetary 
magnetic field data. Munakata et al. (2006) and Kuwabara et al. (2009) obtained an 
improved, more realistic density distribution, derived from the Fokker-Planck equation for an 
isotropic phase-space density. This model used an expanding vacuum cylinder as a magnetic 
cloud model. From the analysis, they deduced the cross-field diffusion coefficient in an 
interplanetary magnetic cloud for high-rigidity cosmic rays observed by muon detectors. 

Recently, the role of gradient- curvature drift for cosmic ray penetration into interplan- 
etary magnetic clouds, which are often observed as magnetic flux ropes, was investigated by 
Krittinatham & Ruffolo (2009). They developed an analytical model of the global magnetic 
field structure of a curved interplanetary magnetic flux rope with its roots connecting to the 
Sun. Cosmic ray particles guiding center orbits were traced in the model. It was shown that 
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gradient-curvature drift plays a significant role for 1 GeV cosmic ray penetration into the 
magnetic flux rope. They emphasized that pitch angle scattering and perpendicular diffusion 
are also important. 

All the models described above only considered guiding center motion and ignored the 
gyration of cosmic ray particles, namely, the Larmor radius of the cosmic ray particle is 
assumed to be infinitely small compared with the flux rope radius. However, in the case of 
high-rigidity cosmic rays inside an interplanetary magnetic flux rope, the assumption is not 
always valid. For example, the rigidity of cosmic rays detected by muon detectors is roughly 
60 GV (e.g., Kuwabara et al. 2004). Assuming a typical magnetic flux rope near the Earth 
with a radius of 0.1 AU and an axial magnetic field intensity of 20 nT (Marubashi 1997), 
the ratio of the Larmor radius to the flux rope radius is then about 0.669, which is not 
much less than unity. This implies that the assumption described above is not valid for the 
high-rigidity cosmic rays inside an interplanetary magnetic flux rope. Therefore, we should 
examine the effect of a finite Larmor radius on cosmic ray penetration into an interplanetary 
magnetic flux rope. 

The purpose of this article is to discuss the effect of a finite Larmor radius on the 
penetration process for cosmic rays into an interplanetary magnetic flux rope from the outside 
by considering the gyration of cosmic ray particles. In Section 12. H we derive equations 
describing the behavior of cosmic ray particles inside a magnetic flux rope without magnetic 
field irregularities on the basis of the Newton-Lorentz equation of each cosmic ray particle. 
The numerical simulations of the cosmic ray penetration into the magnetic flux rope with 
magnetic field irregularities are described in Section [2^21 Sections |3] and H] are devoted to the 
discussion and the conclusions, respectively. 



2. The Model 

2.1. Analytical Model 

According to pitch angle scattering theory, a charged particle is the most effectively 
scattered by a magnetic field irregularity having a spatial scale comparable to the Larmor 
radius of the particle. As discussed in Section [H the Larmor radius of a high-rigidity cosmic 
ray particle, which can be observed using a muon detector, is comparable to the spatial scale 
of a magnetic flux rope near the Earth. This means that the high-rigidity cosmic rays inside 
the flux rope will not be scattered effectively by a small-scale magnetic field irregularity, 
which has a spatial scale much less than that of the flux rope. For this reason, we ignore 
cosmic ray scattering by such an irregularity in the analytic model. This model can exclude 
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an effect of diffusion and consider only effects of gyration and gradient-curvature drift. 

Generally, an interplanetary magnetic flux rope expands as it propagates through in- 
terplanetary space. However, we ignore flux rope expansion because its rate is many orders 
of magnitude less than the speed of a cosmic ray particle. We consider a non-expanding 
cylindrical magnetic flux rope model with radius Rq (e.g., Marubashi 1997). The model is 
defined in cylindrical coordinates {r,ip,z) as 

B = B^e^ + 5,e„ (1) 

B^ = sBMap), (2) 

B, = BoMap), (3) 

P=^ (0<p<l), (4) 

where Bq is the magnetic field intensity along the flux rope axis, and Jq and Ji are zeroth- 
order and first-order Bessel functions of the first kind, respectively; s = 1 and s = — 1 also 
correspond to parallel and antiparallel types of flux rope, respectively. Parallel/antiparallel 
flux ropes are those with electric current flowing parallel/antiparallel to the magnetic field 
in the flux rope. The constant a ~ 2.40483 is the smallest positive number of the zero point 
of the zeroth-order Bessel function of the first kind; that is, Jo{a) = 0. 

The Newton-Lorentz equation for a cosmic ray particle inside a flux rope is written as 

where m, g, v, 7, and c are the mass, electric charge, velocity, Lorentz factor of a cosmic 
ray particle, and the speed of light, respectively. A solution of Equation ([5]) is generally 
described by one parameter t. However, as the solution mathematically defines a curved line 
on the spherical surface in the velocity space, it is possible to choose another parameter as 
long as one-to-one correspondence between the parameter and the time t is confirmed. Since 
a distance r from the flux rope axis to the particle position is described as 

t 
Vrdt + ro, (6) 

to 

where to and ro are constants, r is a monotonous function of time t and there is one-to-one 
correspondence between r and time t as long as the sign of Vr does not change. If the sign of 
Vr changes, we divide the curved line at the time oi Vr = into separate segments. For each 
segment, there is one-to-one correspondence between r and time t because the sign of Vr 
does not change. As a result, r instead of the time t can be applied as a parameter for each 
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segment of the curved line. We can rewrite Equation (|5]) for each segment as an ordinary 
differential equation of r or p (Equation (j4])) by using the relation dr = Vrdt derived from 
Equation (^ as 

du uxb .„. 

dp to 

where u and b are the normalized velocity and magnetic field, respectively, defined as 

u = v/\v\, b = B/Bo. (8) 

/o is defined as 

^ m-f\v\c ^ Rl .g^ 

qBoRo Rq 
where Rl is the Larmor radius of a cosmic ray particle for the magnetic field at the flux rope 
axis. Here, /o corresponds to the ratio of the Larmor radius of a cosmic ray particle to the 
flux rope radius. As the normalized magnetic field in the flux rope b depends only on r or p, 
all segments of the curved line are described by the same equation (Equation ([7])). Also as 
two consecutive segments pass the same point in the phase space at the time of Vr = 0, the 
solutions of Equation ([7]) for all segments should be the same in velocity space. Therefore, in 
the flux rope, the solution of Equation (^^ together with a condition l^l = 1 and the solution 
of Equation (^ depict the same curved line in the velocity space. Equation l^^ shows that 
the behavior of cosmic rays inside the flux rope is regulated only by parameter /q. 

The r, (p and z components of Equation (|7j) are expressed as 

(10) 
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(12) 
The normalized velocity u can be written in spherical coordinates {O^cj)) as 

u^ = sin 6^ cos 0, u<^ = sin 6' sin 0, u^ = cos6', (13) 

where 9 and are measured from the z axis and r axis, respectively (see Figured]). From 
the definition of the flux rope, the normalized magnetic field components b^p and b^ are 



and 
where 

and 
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b^ = sJi{ap), b-, = Jo{ap), (14) 

respectively. 

We set the boundary conditions as Ur = sin6'ocos0o; u^ = sin^osin^o, and u^ = cos 6*0 
at p = po < 1- Equations flTO|) . f lTTj) . and ( !T2|) can be solved analytically. The solutions are 
expressed as 

ul = 1 ^ (po sin 6o sin 0o - C^) - (cos 610 - C^) , (15) 

Po sin 6^0 sin 00 - C'vp ^,r\ 

u^ = ^, 16 

P 

u^ = cos9o-C^, (17) 

^ ^ pJi(ap) -poJi(apo) ,^g, 

'^ afo 

fj ^ s [Map) - MaPo)] QgN 

These solutions include gyration and gradient-curvature drift effects. Equation f lT^ shows 
that the particle position p should satisfy the relation: 

(po sin 60 sin 00 - C^? + p\cos 60 - C,f - p' < 0. (20) 

This means that only particles having the boundary condition ^o and 0o at po can reach the 
region where Equation (!20|) is satisfied. 

To examine the validity of the transformation from Equation (jS]) to Equation (jTj), we 
compare the analytical solutions with a numerical simulation, in which a cosmic ray particle 
is directly traced in the magnetic flux rope. We integrate Equation ([5]) and v = dr/dt by 
using the Buneman-Boris Method. Needless to say, the numerical simulation includes both 
gyration and gradient-curvature drift effects. Hereafter, we consider a parallel-type magnetic 
flux rope model (i.e., s = 1) without loss of generality. We use Rq = 0.1 AU and i?o = 20 nT 
as typical values near the Earth (Marubashi 1997). Figure [2] shows the particle trajectories 
on the r — (f plane. A dashed circle with a radius of 1 in the figure is the edge of flux rope. 
The solid curves starting from the point {x,y) = (1,0) and {x,y) = (—1,0) are the particle 
trajectories of 60 GV cosmic rays, which correspond to /o = 0.669, and those of 6 GV cosmic 
rays, which correspond to /o = 0.0669, respectively, for various boundary conditions 9q and 
00 at Po = 1. The starting points are not important for the results because the magnetic 
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field in the flux rope is symmetrical with regard to the z axis. All these particles entering 
from the outside of the flux rope are escaping from the flux rope within one gyro motion. A 
thick spiral curve in the figure shows a trajectory of a 6 GV trapped cosmic ray particle for 
randomly selected boundary condition d^ and 00 at po = 0.5. The trapped particle trajectory 
obviously shows a gradient-curvature drift. The drift velocity components are in the y? and 
z directions as calculated from the definition of the flux rope. 

Figure [3] shows the three components of normalized velocity u against the normalized 
distance p from the flux rope axis for 60 GV cosmic rays. Panels a), b), and c) show the r, 
(y9, and z components, respectively. Gray plus signs show results of numerical simulation for 
various penetration angles ^o and c^o at po = 1, corresponding to those in Figure [2l Solid 
lines show the analytical solutions for the same parameters used in the numerical simulation. 
It is evident that the analytical solutions exactly trace the simulation results. It is noticed 
that unphysical solutions, in which Ur is an imaginary number, are also drawn by dashed 
lines in Panels b) and c) in the figure. 

Figure m also shows the three components of normalized velocity against the normalized 
distance for a 6 GV trapped cosmic ray particle. Gray plus signs show results of numerical 
simulation, represented by the spiral curve in Figure [2j Solid lines are the analytical solutions 
for the same parameters used in numerical simulation. It is also evident that the analytical 
solutions exactly trace the simulation results. It is again noticed that unphysical solutions, 
in which u^ is an imaginary number, are also drawn in Panels b) and c) in the figure. The 
trapped cosmic ray particle obviously gyrates and drifts in the magnetic fields, as described 
above. The trajectories of these particles in velocity space show oscillatory motions on the 
lines described by our analytical solutions. Figures [21 [3l and H] indicate that our analytical 
solutions represent exact particle motion including particle gyration and gradient- curvature 
drift in the magnetic flux rope. 

Hereafter, we consider cosmic rays penetrating a flux rope from the outside; that is, ^o 
and 00 at po = 1 as the boundary condition. Not all particles reach the axis p = in Figure 
[3] because only particles with specific angles 6^0 and 0o satisfy Equation fl20l) at p = 0. A 
cosmic ray particle entering the flux rope at angles Qq and 00 moves along the lines shown 
in the Panels b) and c) in Figure [3] until it reaches the innermost point where the left-hand 
side of Equation (!20l) is zero. The particle then returns along the same line to the edge of 
the flux rope (p = 1). The (/? and z components of the velocity of the escaping particle are 
the same as those of the entering particle if p is the same, and the r component of velocity 
is in the opposite direction and has the same magnitude when the particle is entering and 
escaping. Panel c) in Figure [3] and Equation flTTl) show that the z component of the velocity 
decreases when the particle is entering and increases when it is escaping. This is because 



the z component of the Lorentz force is always negative (positive) when cosmic ray particles 
move from (to) the edge of the flux rope to (from) the innermost point. We find from 
Equation ( IT^ that the slopes of the trajectories in the panel c) in Figure E] are steeper for 
lower-rigidity cosmic rays and less steep for higher-rigidity cosmic rays (not shown in the 
figure) . 

Cosmic ray particles penetrating from the outside of a flux rope must have m^ > at 
the edge of a flux rope, po = 1- From this fact and Equations (IT^ to ( 1T^ . we find that the 
cosmic ray particles penetrating from the outside of the flux rope cannot exist in the region 
that satisfies 

l-(psin^sin0 + C^)^-(cos^ + a)^ < 0, (21) 

when we call it a trap region. Figure O shows the trap region for /q = 0.0669. Panels a) 
and b) are the trap region for p = 0.85 and 0.95, respectively. These panels show that the 
solid angle of the trap region depends on p. The solid and dashed lines in Figure [6] show p 
dependence of the solid angle of the trap region normalized by 4tt steradian for /o = 0.0669 
and /o = 0.669, respectively. The solid angle of the trap region is larger (smaller) for 
smaller (larger) p for a given /q, and is Air steradian in the inner region of the flux rope for 
/o = 0.0669. If cosmic ray particles are in the trap region, the particles cannot escape from 
the flux rope and are trapped. 

We calculate the forbidden region to which a cosmic ray particle cannot arrive from the 
outside of the flux rope for any 9q and 0o using Equation ( 12U]) . Figure [7] depicts the forbidden 
region. The abscissa and ordinate are /o and p, respectively. The figure shows that low- 
rigidity cosmic rays, which can be observed by a neutron monitor, usually cannot arrive at 
the inner region of the flux rope. For example, in the case of /o = 0.0669, which corresponds 
to a 6 GV cosmic ray particle, the particle cannot penetrate the region for p < 0.769. Thus, 
the low-rigidity cosmic rays cannot reach the inner region of a magnetic flux rope by only 
gyration and gradient- curvature drift. 

However, many events of cosmic ray decreases observed by neutron monitors suggest 
that cosmic rays exist in the inner regions of flux ropes. Therefore, other mechanisms for 
penetration into the inner region of flux rope by low-rigidity cosmic ray particles are required. 
For these particles, /o is much less than unity and they can be scattered effectively by small- 
scale magnetic field irregularities in the flux rope. For this reason, we consider small-scale 
magnetic field irregularities in the interplanetary magnetic flux ropes in Section 12.21 

On the other hand, since the Larmor radius of high-rigidity cosmic ray particles is 
comparable to the flux rope radius (e.g., /o = 0.669 for 60 GV), these particles can reach the 
inner region of the flux rope from the outside by gyration. While the high-rigidity cosmic 
rays can reach the inner region of the magnetic flux rope by only gyration, the effect of the 
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small-scale magnetic field irregularities on the cosmic ray penetration into the flux rope is not 
clear. Therefore, we also investigate the effect of the small-scale magnetic field irregularities 
in the interplanetary magnetic flux ropes for high-rigidity cosmic rays in Section 12.21 

In this study, we only considered the case of a flux rope located around the Earth. Since 
a magnetic flux rope is ejected from the Sun and propagates through interplanetary space 
toward the Earth while expanding, we must consider the flux rope near the Sun. According 
to the expanding cylindrical force-free magnetic flux rope model (Shimazu & Vandas 2002; 
Berdichevsky et al. 2003), the flux rope radius Rq is directly proportional to the time t 
that has passed since the flux rope ejection at the Sun, and the magnetic field intensity on 
the flux rope axis Bq is inversely proportional to time t squared. From these dependences 
of Rq and Bq, and Equation (^, the parameter /q is directly proportional to time t. This 
means that the parameter /q near the Sun is much less than that around the Earth; for 
example, /o = 0.0669 for a typical flux rope located 0.1 AU from the Sun and 60 GV cosmic 
ray particles. This is the same situation as having low-rigidity cosmic ray particles around 
the Earth because the behavior of cosmic rays inside a flux rope is regulated by only the 
parameter /q. Therefore near the Sun, even high- rigidity cosmic rays cannot penetrate the 
interior of the flux rope by gyration alone, and the role of magnetic field irregularities must 
also be considered. 



2.2. Numerical Model 

To examine the role of magnetic field irregularities in the penetration of cosmic rays into 
interplanetary magnetic flux ropes, we perform numerical simulations in which cosmic ray 
particles are directly traced in the turbulent magnetic field. The direct simulations solving 
the Newton-Lorentz equation in the turbulent magnetic field were recently performed by 
some authors, for instance, in studies of cross-field diffusion (e.g., Michalek & Ostrowski 
1997; Giacalone & Jokipii 1999; Mace et al. 2000; Qin et al. 2002; Zimbardo et al. 2006), 
pitch angle diffusion (Qin & Shalchi 2009), suppression of particle drift (Minnie et al. 2007), 
and solar energetic particle transport (Qin et al. 2004). 

We trace the cosmic ray particles in the magnetic flux rope with the addition of magnetic 
field irregularities. Magnetic field irregularities are defined as the superposition of many 
sinusoidal waves, 

Btotai = B + 6B, (22) 

TV 

6B = Y, CnBCiK) exp [t (fc„ ■ r + z^„)] , (23) 
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C{knf= — F{kn)Akn y^F{kn)Akn , (24) 

(25) 

where B = \B\, 6B = \SB\, and kn = |fc„|. Vector ^„ is perpendicular to fc„ to assure 
divBtotai = 0. Equation ( 125|) defines a Kolmogorov-type wave number spectrum of magnetic 
field irregularities. These definitions are similar to those of previous studies (e.g., Giacalone & 
Jokipii 1999). The minimum and maximum wave number used in the simulation correspond 
to wave length of 2Ro (diameter of flux rope) and 2 x 10~^i?o, respectively. Ic is a correlation 
length of the turbulence spectrum and is set as 0.5i?o- One sinusoidal wave is determined by 
four random numbers; two for a direction of wave number vector, one for a phase of wave, and 
the other for a polarization. The 255 sinusoidal waves are superimposed on the background 
magnetic field. We confirmed that these waves are sufficient to reproduce a diffusion process 
and that results are insensitive to the number of waves. We use a value [SB/B]"^ = 0.01 as 
an intensity of magnetic field irregularities in this study, which is consistent with values of a 
small percent found in some magnetic clouds (e.g., Subramanian et al. 2009). 

In the simulation, the background magnetic field B is the fiux rope magnetic field 
defined as Equations ([I]) to (jl]). We consider a parallel- type magnetic fiux rope model with 
the typical values of -Rq = 0.1 AU and Bq = 20 nT near the Earth. For a time dependence of 
fiux rope radius and magnetic field, we use the model by Shimazu & Vandas (2002), which 
can very effectively reproduce observations. In this model, when the flux rope is located 
near the Sun, for instance, at 0.1 AU from the Sun, Rq and Bq will be 0.01 AU and 2,000 
nT, respectively. 

Cosmic ray particles are traced by the Buneman-Boris Method. The number of particles 
M is 10,000. All these particles are initially located at the flux rope edge with randomly 
distributed penetration angles. When a particle escapes from the edge, a new particle is set 
up immediately at the edge and we trace its orbit. A different set of the random numbers 
for magnetic fleld irregularities are selected for each particle. 

Figure [8] depicts the simulated particles on the r — ip plane for a case where the flux 
rope is located at 0.1 AU from the Sun and the cosmic ray rigidity is 60 GV. In this case, the 
parameter /q is 0.0669. Solid and dashed circles in the flgure show the flux rope edge and the 
boundary of the forbidden region shown in Figure [TJ respectively. The gray tiny points are 
cosmic ray particles. Panels a), b), c), and d) are snapshots for the time t = 2.5, 250, 750, 
and 2,500 ou~^, where Up is the proton cyclotron frequency {up = 0.479 sec~^ in this study). 
Panel a) shows that almost all cosmic ray particles are located outside of the forbidden region 
and the effect of magnetic fleld irregularities is small. However, after approximately nine 
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minutes later, in Panel b), the effect of magnetic field irregularities clearly appears, and a 
large number of particles penetrate into the forbidden region in the flux rope. In Panel c), 
some particles reach the center of the flux rope. In Panel d), cosmic ray particle distribution 
is nearly uniform. Penetration into the forbidden region is due not to drift motion but to the 
diffusion process caused by particle scattering and magnetic field line random walk (Jokipii 
& Parker 1968). 

Figure [H] shows a time evolution of a mean radial distance of 60 GV cosmic ray particles 
in the typical flux rope, whose parameters are Rq = 0.1 AU and Bq = 20 nT near the Earth. 
The mean radial distance at time t is calculated as 

^ k M 

^ ^ j=0 i=l 

where k is simulation time step number until time t and pij is a distance from the flux rope 
axis of the zth cosmic ray particle at the simulation time step j. In this definition, the 
entire simulation time step up to time t is used to calculate (p). This definition is based 
on a continuous cosmic ray penetration from the outside of the flux rope. The thin dotted 
horizontal line in Figure |9] shows the result if the particles are uniformly distributed in the 
flux rope, (p) = 2/3. Thick and thin dashed lines show the cases without magnetic field 
irregularities near the Sun (/o = 0.0669) and Earth (/o = 0.669), respectively. Thick and 
thin solid lines show the cases with magnetic field irregularities for near the Sun and Earth, 
respectively. We find from the figure that the mean radial distance immediately reaches 
steady state values in all cases. In the case near the Sun (thick solid line), the mean radial 
distance does obviously decrease and is asymptotic to the uniform distribution line as the 
time elapses. This means that the cosmic ray particles penetrate into the inner region of the 
flux rope by diffusion, as shown in Figure El A mean radial distance for the case near the 
Earth (thin lines) is already small at the beginning of the calculation compared to that for 
the case near the Sun. This is because cosmic ray particles are able to reach the inner region 
of flux rope by gyration because of the large /o parameter. However, the decrease of the 
mean radial distance is small even if the time elapses. We find from the figure that it takes 
more time to penetrate into the flux rope by diffusion; thus the gyration is more effective 
than the diffusion for high /q such as near the Earth. The diffusion is effective when the flux 
rope is located near the Sun (small /o). 

Figure ITU] shows the mean radial distance for 6, 18, 24, and 48 GV cosmic ray particles 
for the flux rope located near the Earth. The simulation run time is 10 times longer than 
that in the previous case shown in Figure |9l The solid, dashed, dash-dotted, and long-dashed 
lines are the results for 6, 18, 24, and 48 GV cosmic ray particles, which correspond to /o = 
0.0669, 0.201, 0.268, and 0.535, respectively. The dotted horizontal line shows the mean 
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radial distance for which the particles are uniformly distributed in the flux rope. The mean 
radial distances for all cases clearly decrease with time. We find from comparison between 
the thick solid line in Figure M and the solid line in Figure [TD] that cosmic ray particles 
attain the same distribution in the flux rope if the parameter /q is the same even when the 
magnetic field irregularities exist. The time taken to attain this distribution is proportional 
to the flux rope radius. This fact means that some distribution of the cosmic ray density 
can be attained in a short time when the flux rope is located near the Sun. It is interesting 
that the mean radial distances drop under the uniform distribution line in all cases except 
for the 6 GV case. This means that a density in the inner region of a flux rope becomes 
larger than the density in outer region. This result shows the significance of the effect of 
finite Larmor radius because this kind of density distribution cannot be attained by only the 
diffusion process. 

Figures [11] to [T3] show the calculated cosmic ray radial density distributions inside the 
flux rope located near the Earth for the cases of 6, 18, 24, and 48 GV cosmic ray particles, 
respectively. All the simulation time steps up to time t are used to calculate the density 
distributions at time t. The dashed lines in the figures are density distributions in the 
steady state, which are calculated by simulations without magnetic field irregularities. Ten 
solid lines from the bottom to the top are density distributions calculated by simulations 
with magnetic field irregularities for time t = 2,500 to 25,000 u~^ with an interval of 2,500 
o;"^, respectively. Dotted vertical lines in Figures ITT] and [T2] show the forbidden distance, 
which are borders of the forbidden region calculated theoretically in Section 12.11 and shown 
in Figure [3 

We can see from Figure [TTl which shows the case for 6 GV cosmic ray particles, that 
the density distribution inside and outside the forbidden region is greatly different. Out of 
the region, there is a density bump around p = 0.8 in the later time. The density out of 
the forbidden region increases drastically with time and the bump is formed just outside the 
region. The bump is formed in the following processes. The particles penetrating from the 
outside of the flux rope edge are scattered by small-scale magnetic field irregularities between 
the forbidden region and the flux rope edge, and changed their pitch angles. These pitch 
angle changes can make the particles fall into the trap region, which is expressed as Equation 
f l?I]) and drawn in Figure El As a result, the particles are trapped between the forbidden 
region and the flux rope edge. As the solid angle of the trap region near the forbidden region 
is larger than that near the flux rope edge (see Figure [6]), cosmic ray particles tend to fall in 
the trap region near the forbidden region, not near the edge. Accumulating these trapped 
particles, the bump in the density distribution is formed. Moreover, there is a small density 
drop just inside the flux rope edge. The drop is formed in the following processes. Since the 
solid angle of trap region is very small just inside the flux rope edge as shown in Figure |6l 
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nearly particles cannot be trapped in this region. Therefore, the cosmic ray particles in this 
region are distributed like the steady state distribution without magnetic field irregularities 
(dashed line in Figure [TT]) . As the density bump is formed near the forbidden region, this 
fills in much of the density drop just inside the edge of the flux rope, which is why the 
small density drop is formed just inside the flux rope edge. A major role of magnetic fleld 
irregularities outside of the forbidden region is to make the particles fall into the trap region 
by pitch angle scattering. Outside of the forbidden region, the effect of a flnite Larmor radius 
is significant for density distribution. Inside of the forbidden region, the trapped particles 
just out of the forbidden region diffuse into the forbidden region by means of spatial diffusion 
due to scattering by magnetic field irregularities and a magnetic field line random walk, which 
is why the density distribution in the forbidden region decreases monotonically toward the 
fiux rope center. Thus, a major role of magnetic field irregularities inside the forbidden 
region is spatial diffusion toward the flux rope center. 

Figure [12] shows the results for the case of 18 GV cosmic rays particles. The bump in 
the density distribution appears while the peak of the bump is located closer to the flux 
rope axis compared to the case of 6 GV cosmic ray particles. This is because the effect of 
the Larmor radius for 18 GV cosmic ray particles is greater than that for 6 GV. A density 
distribution in the forbidden region can be interpreted as the diffusion process as in the case 
of the 6 GV cosmic ray particles mentioned above. 

The density distributions for 24 GV and 48 GV cosmic ray particles are shown in 
Figures [12] and [TU, respectively. These medium- to high-rigidity particles have no forbidden 
regions. The density distribution for 24 GV cosmic ray particles is almost flat in the inner 
region although the density peak is formed deep inside in the later time. Although no bump 
distribution is found in the 48 GV case, the density in the inner region is larger than that in 
the outer region. Since density in the inner region is smaller than that in the outer region 
for the density distribution derived from only the diffusion process (Cane et al. 1995), the 
distributions derived from our simulations, which include an effect of a finite Larmor radius, 
are reverse distributions compared to those derived from the diffusion process only. In these 
rigidity ranges, a major role of magnetic field irregularities is to make the particles fall into 
the trap region by pitch angle scattering. These facts show that the effect of a finite Larmor 
radius is significant, and we must consider that effect. 



3. Discussion 

In this study, we used an isotropic magnetic field irregularity model in an interplanetary 
magnetic fiux rope. However, the irregularities in the magnetic fiux rope are likely to be 
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anisotropic (Leamon et al. 1998; Narock & Lepping 2007). Anisotropic magnetic field ir- 
regularities cause somewhat large perpendicular diffusion compared to that due to isotropic 
magnetic field irregularities (Giacalone & Jokipii 1999). The intensity of magnetic field irreg- 
ularities {[5B / B]'^ = 0.01 in this study) may also be a little different. Therefore, our results 
may be modified somewhat quantitatively by the anisotropic magnetic field irregularities and 
its intensity. In our model, we also used a minimum and maximum wave number, and the 
correlation lengths of the turbulent spectrum were set as the specified values. These values 
may be a little different from the real situation. We performed many simulations for other 
parameters and confirmed that they did not cause the results to change qualitatively. 

Generally, the magnitude of Forbush decreases is highly variable between several to a 
couple of tens percent for neutron monitor observations (Cane et al. 1996). However, our 
model shows a much larger decrease (see Figure [TT]). The reason for this difference is that 
simulation time is shorter than propagation time from the Sun to the Earth. The simulation 
time in our model is 25,000 oo"^, corresponding to 0.6 days, because our simulation is time- 
consuming. In reality, however, the time when a flux rope propagates from the Sun to the 
Earth is about 1.5 to 2.5 days. If the simulation runs for a longer time, the cosmic ray 
density in the flux rope will increase and magnitude of decrease will be smaller. We have 
now investigated more realistic models with longer simulation time. 

We should discuss the effects of the flux rope curvature and magnetic field configuration 
for a realistic situation. First, we discuss the effect of curvature of the axis of an interplan- 
etary magnetic flux rope. In the curved flux rope, a cosmic ray particle can have a radial 
component of gradient-curvature drift velocity and a cosmic ray can penetrate into the flux 
rope, as discussed by Krittinatham & Ruffolo (2009). They demonstrated that cosmic ray 
particles can penetrate into the outer part of a flux rope by means of drift and penetrate 
into the inner part of a flux rope by diffusion. This is similar to our result, which is the pen- 
etration by diffusion into the inner part of the flux rope. On the other hand, the mechanism 
of penetration into the outer part is different from theirs. This difference is caused by an 
effect of the finite Larmor radius. Our result shows that cosmic ray particles can penetrate 
into the outer part of a flux rope via the effect of the finite Larmor radius. Cosmic ray 
penetration into the curved flux rope with this effect should be further investigated and will 
advance. Second, we discuss the magnetic field configuration. The roots of magnetic field 
lines of some interplanetary magnetic flux ropes may be detached from the solar surface and 
connected to interplanetary magnetic fields (Gosling et al. 1995; Cane et al. 2001). In this 
situation, cosmic ray particles for a low /q parameter may be able to penetrate into flux 
ropes along interplanetary magnetic field lines because the Larmor radius is small compared 
with the flux rope radius. This may also fill the inner region with cosmic ray particles. 
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An interplanetary magnetic flux rope is often accompanied by a shock wave in front 
of it. This shock wave can accelerate particles and be one of an origin of solar energetic 
particles. These energetic particles can also penetrate across the flux rope boundary. The 
penetration of these particles is, mainly, diffusion and gradient-curvature drift, and the effect 
of the finite Larmor radius is small because the particles' energy is several MeV to GeV and 
the Larmor radius is negligibly small compared to the flux rope radius. 



4. Conclusions 

We studied processes for the penetration of cosmic rays into an interplanetary magnetic 
flux rope, focused on an effect of the finite Larmor radius and magnetic field irregularities. 
We considered a cylindrical magnetic flux rope as a model of an interplanetary magnetic 
flux rope. First, we derived an analytical solution for cosmic ray behavior in a magnetic flux 
rope under the assumption there is no scattering by small-scale magnetic field irregularities. 
The solution showed that cosmic ray behavior is determined by only the parameter /q; 
that is, the ratio of the cosmic ray Larmor radius at the fiux rope axis to the fiux rope 
radius. We showed that cosmic rays cannot penetrate into the inner region of a fiux rope 
by only gyration and gradient-curvature drift in the case of small /q. Next, we performed 
a numerical simulation of cosmic ray penetration into an interplanetary magnetic flux rope. 
The simulation included an effect of small-scale magnetic field irregularities. The results 
showed that cosmic rays can penetrate into the magnetic flux rope even in the case of small 
/o by the effect of small-scale magnetic field irregularities. The result deduced from the 
simulation is that the high-rigidity cosmic ray penetration is dominated by diffusion when 
the flux rope is located near the Sun while dominated by gyration when it is located near 
the Earth. This simulation also showed that a cosmic ray density distribution for the case 
of moderate to large /o parameter is greatly different from the distribution deduced from a 
guiding center approximation shown in previous studies. Our results show that the effects of 
the finite Larmor radius and magnetic field irregularities cause the bump or reverse cosmic 
ray density distribution inside the interplanetary magnetic flux rope. 

The conclusion of this study is that an effect of the finite Larmor radius is significant for 
the cosmic ray penetration into an interplanetary magnetic flux rope if the Larmor radius of 
cosmic ray particles in the flux rope is not negligible compared to the flux rope radius. 

We would like to thank the anonymous referee for the useful comments on improving 
the manuscript. 
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Fig. 1. — Definition of tlie coordinate system for the cosmic ray velocity v in the interplane- 
tary magnetic flux rope. The origin is located on the cosmic ray particle. Here, 6 is measured 
from the z axis, which is parallel to the flux rope axis, while is measured from the r axis, 
which leads from the center of the flux rope to the particle. The spiral line represents a 
magnetic field line of the flux rope. 
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Fig. 2. — Particle trajectories on the r — ip plane observed in the numerical simulation. The 
dashed circle with a radius of 1 is the edge of the flux rope. The solid curves starting from 
the point {x,y) = (1,0) and {x,y) = (—1,0) are the particle trajectories of 60 GV cosmic 
rays and 6 GV cosmic rays for various boundary conditions 6q and (pQ at Pq = 1, respectively. 
The thick spiral curve shows the trajectory of a 6 GV trapped cosmic ray particle (see text) 
for randomly selected boundary condition 6q and 0o at po = 0.5. A 0.1 AU radius flux rope 
with a 20 nT axial magnetic fleld is used for typical values near the Earth. 
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Fig. 3. — Three components of the normahzed velocity for 60 GV cosmic ray particles for 
various penetration angles 6*0 and 0o at po = 1, which correspond to those in Figure |2j The 
components are depicted against the normalized distance from the flux rope axis. The edge 
of the flux rope is at p = 1. Panels a), b), and c) show the r, ip, and z components of 
the velocity, respectively. The solid lines are expressed as Equations (IT5l) . ( IT6l) . and ( IT71) . 
Unphysical solutions, in which Ur is an imaginary number, are also drawn by dashed lines. 
The gray plus signs show results of the numerical simulation. 



-21 - 



a) u^ 



b) u 



0.5 







-0.5 



■1 



n 1 r 






C) U^ 




0.3 0.4 0.5 0.6 0.3 0.4 0.5 0.6 0.3 0.4 0.5 0.6 



Fig. 4. — Three components of the normahzed velocity for a 6 GV trapped cosmic ray particle 
for angle Oq and (f)Q at po = 0.5, which correspond to those in Figure O The components 
are depicted against the normalized distance from the flux rope axis. The solid and dashed 
lines, and gray plus signs, are the same as those in Figure [31 
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Fig. 5. — Trap region (gray) for /o = 0.0669. Panels a) and b) are the trap region for 
p = 0.85 and 0.95, respectively. The solid angle of the trap region depends on p. 
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Fig. 6. — p dependence of the solid angle of the trap region normalized by 4tt steradian. 
Solid and dashed lines are for /o = 0.0669 and /o = 0.669, respectively. The solid angle of 
the trap region is larger for smaller p for a given /o, and is Att steradian in the inner region 
of the flux rope for /q = 0.0669. In this region, all particles are trapped. 
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Fig. 7. — Forbidden region (gray) where cosmic rays penetrating from the outside cannot 
enter. The abscissa and ordinate are, respectively, the parameter /q and distance from the 
flux rope axis normalized by the flux rope radius. 
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Fig. 8. — Cosmic ray particle distribution on the r — ip plane for the case that the flux rope is 
located at 0.1 AU from the Sun and that the cosmic ray rigidity is 60 GV. Solid and dashed 
circles show the flux rope edge and the boundary of the forbidden region, respectively. Panels 

a), b), c), and d) are for the time t = 2.5, 250, 750, and 2,500 (oUp^). 
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Fig. 9. — A mean radial distance calculated from the simulation results. Thick and thin 
lines are the results for the case near the Sun and Earth, respectively. Solid and dashed 
lines are the results with and without magnetic field irregularities, respectively. The dotted 
horizontal line is a mean radial distance for the uniformly distributed case. 
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Fig. 10. — A mean radial distance calculated from the simulation results for the case of near 
the Earth. Solid, dashed, dash-dotted, and long-dashed lines are the results for 6, 18, 24, 
and 48 GV cosmic rays, respectively. Flux rope parameters are i?o = 0.1 AU and Bq = 20 
nT. The dotted horizontal line is the same as in Figure |9l 
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Fig. 11. — Cosmic ray density distributions for 6 GV cosmic rays calculated by simulations. 
Flux rope parameters are the same as the ones used in Figure [TOl The dashed line is the 
distribution for the case without magnetic field irregularities. Ten solid lines from the bottom 
to the top show the time evolution of the density distribution for the case with magnetic 
field irregularities. The region inside the dotted vertical line shows the forbidden region 
calculated from the analytical model. The density n is normalized to unity at the edge of 
the flux rope. 
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Fig. 12. — Same as Figure ITTl except the cosmic ray rigidity is 18 GV. 
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Fig. 13. — Same as Figure [TH except the cosmic ray rigidity is 24 GV. There is no forbidden 
region for this rigidity cosmic ray particle. 
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Fig. 14. — Same as Figure [TH except the cosmic ray rigidity is 48 GV. There is no forbidden 
region for this rigidity cosmic ray particle. 



